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ABSTRACT 

We present a model for determining the mass ratio in interacting binaries by 
directly fitting the observed spectrum with synthetic spectra. We make direct use 
of NextGen model atmospheres intensities which are the most comprehensive and 
detailed models available for cool stars. We fully take into account the varying tem¬ 
perature and gravity across the secondary star’s photosphere, by incorporating the 
synthetic spectra into the secondary star’s Roche geometry. As a result, we determine 
the exact rotationally broadened spectrum of the secondary star and so eliminate the 
need for a limb-darkening law, and the uncertainties associated with it. 

As an example we determine the mass ratio for the well studied soft X-ray tran¬ 
sient Nova Sco 1994. In order to obtain a more accurate determination of the mass 
ratio, which does not depend on assumptions about the rotation profile and limb- 
darkening coefficients, we use our model to compute the exact rotationally broadened 
model spectrum, which we compare directly with the observed intermediate resolu¬ 
tion spectrum of Nova Sco 1994. We determine the mass ratio of Nova Sco 1994 to 
be 0.419±0.028 (90 percent confidence), which is the most accurate determination 
of the binary mass ratio in an X-ray binary. This result combined with the binary 
mass function and inclination angle gives a refined black hole mass of 5.99±O.42M0 
(90 percent confidence). We also perform simulations which show that, for an F-type 
secondary star, the standard rotation profile with zero and continuum value for the 
line limb-darkening coefficient gives a value for q that brackets the value found using 
the full geometrical treatment 

Key words: binaries: - close - stars: fundamental parameters, individual: X-Ray 
Nova Sco 1994 (GRO J1655-40) - X-rays: stars 


1 INTRODUCTION 


It is clear that if one wants to determine accurate binary 
masses for the SXTs, the uncertainties in measnring the 
inclination angle and mass ratio have to be reduced. This 
paper concentrates on reducing the inherent uncertainties 
associated in determining the binary mass ratio. 


The determination of binary masses in interacting bina¬ 
ries where only one component is observable, such as in 
cataclysmic variables (CVs) or low-mass X-ray binaries 
(LMXBs), requires a radial velocity study of the mass loos¬ 
ing star, a knowledge of the binary inclination angle and 
mass ratio. The binary inclination in non-eclipsing binaries 
is usually obtained by measuring the ellipsoidal modulation 
of the late-type star (e.g. see Wilson & Devinney 1971 and 
Shahbaz, Naylor & Charles 1993) and the binary mass ratio 
is normally obtained by measuring the rotational velocity of 
the secondary star (Wade & Horne 1988). This procedure 
has been used for many years to determine binary masses, 
in particular in the subclass of the LMXBs called the soft 
X-ray transients (SXTs), which are characterised by their 
episodic outbursts and contain some of the best stellar mass 
black hole candidates (van Paradijs & McClintock 1995). 


Shahbaz (1998) showed that in principle, one could use 
the shape of the secondary star’s absorption lines to deter¬ 
mine the binary system parameters, such as the inclination 
angle. However, this requires high resolution echelle spec¬ 
troscopy, which for most of the faint SXTs is very difficult, 
even with 10-m class telescopes. Due to the faintness of the 
secondary star in the SXTs, the secondary star’s rotational 
broadening is usually determined by using intermediate res¬ 
olution (^lA) spectroscopy. The secondary stars typically 
have Kot sin i values of 30-100 km/s. With intermediate 
resolution spectroscopy the information about the shape of 
the absorption lines is lost, the only information that can 
be extracted is the amount by which the star’s spectrum is 
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broadened. The procedure commonly used to measure the 
secondary star rotational broadening is to compare it with 
the spectrum of a slowly rotating template star, observed 
with the same instrumental conhguration which has been 
convolved with a limb-darkened standard rotation profile 
(Gray 1992; Collins & Truax 1995 and references within). 
For this one needs to assume a limb-darkening coefficient 
for each spectral line. Usually one either adopts 0, i.e. no 
limb darkening, or the continuum value (which depends on 
the wavelength and the star’s effective temperature). The 
width of the standard rotation profile is varied until an op¬ 
timum match is found with the width of the target spectrum 
(e.g. see Marsh, Robinson & Wood 1994). 

It should be noted that there are many assumptions 
inherent in using this method, primarily due to the use of 
the standard rotation profile. The standard rotation profile 
assumes that the star is spherical and that the line profile 
has the same shape over the entire star. It also assumes a 
wavelength dependent limb darkening coefficient and a ho¬ 
mogeneous stellar surface. Firstly, it should be noted that 
the secondary star in a CV or LMXB substantially fills its 
Roche lobe so is far from being spherical and will thus have 
distorted line profiles the extent of which depends on the 
exact Roche geometry (Shahbaz 1998). Also, the spectrum 
of the secondary cannot be described by a single-star spec¬ 
trum, since both temperature and gravity vary over the pho¬ 
tosphere due to the shape of the star’s Roche-lobe. Further¬ 
more, stellar absorption lines have core limb-darkening co¬ 
efficients much less than the appropriate continuum value 
because the line flux arises from higher regions in the stellar 
atmosphere than the continuum flux; the precise value for 
the line limb-darkening coefficient depends on many param¬ 
eters such as the incidence angle and so requires detailed 
calculations (Collins & Truax 1995). Using zero or the con¬ 
tinuum value introduces a systematic uncertainty of about 
14 percent in the determination of Kot sin i (Welsh et ah, 
1995, Shahbaz 1998) which contributes significantly to the 
accuracy which one can determine the mass ratio and hence 
the binary masses. 

In an attempt to reduce the uncertainties in the binary 
mass ratio derived from intermediate spectroscopy using this 
procedure, we determine the exact rotationally broadened 
spectrum from the secondary star in an interacting binary. 
Our calculations take into account the shape of the sec¬ 
ondary star’s Roche-lobe and the varying temperature and 
gravity of different photospheric segments, which are char¬ 
acterised by individual model atmospheres. By directly com¬ 
paring the model spectrum with the observed intermediate 
resolution spectra of an SXT, we show that we can accu¬ 
rately determine the binary mass ratio. As an example, we 
determine the binary mass ratio for the well studied soft 
X-ray transient. Nova Sco 1994 (=GRO J1655-40). 


2 DETERMINING THE SPECTRUM IN 
INTERACTING BINARIES 

Although there are many codes that use model atmospheres 
to determine the light curves of binary systems (Wilson & 
Devinney 1971, Tjemkes, van Paradijs & Zuiderwijk 1986, 
Shahbaz et ah, 1993, Orosz & Hauschildt 2000) there are 
not many that compute the observed spectrum. Linnell & 


Hubeny (1994, 1996) have a code that computes the syn¬ 
thetic spectra and light curves for binary stars. However, 
since they use Kurucz (1979) model atmospheres to syn¬ 
thesize the spectrum, their code is really only applicable 
for systems with relatively hot stellar components i.e. early- 
type stars. In LMXBs, the low-mass companions are typi¬ 
cally late-type K stars with effective temperatures less than 
5000 K. Also, since the temperature varies across the star’s 
surface due to gravity darkening, there are elements on the 
star’s surface much cooler than the star’s effective tempera¬ 
ture, especially near the inner Lagrangian point (Li; Shah¬ 
baz 1998). Therefore, in order to compute the observed spec¬ 
trum of the late-type secondary star in an LMXB accurately 
and to make the code more general for a range of secondary 
stars, we have written a code that is similar to the code used 
in Shahbaz (1998) but with a full treatment of limb dark¬ 
ening through the addition of synthetic NextGen model 
atmospheres. 


3 THE MODEL 

The equations that determine the basic geometry of an in¬ 
teracting binary system have been known for a long time and 
are well understood. To determine the light curve or spec¬ 
trum from the secondary star in a close binary, one would 
compute the specific intensity for each element of area on 
the star’s surface by solving the radiation transfer equation 
for a given model photosphere. By doing this one eliminates 
the need to use a limb-darkening law which parameterizes 
the specific intensity for different emergent angles, given the 
monochromatic intensity. However, since the computation 
of stellar atmospheres is computationally intensive, approx¬ 
imations are required to determine the expected light curve 
or spectra from an interacting binary system. 

The main geometrical part of the model we use is the 
same as that described in Shahbaz (1998) and is similar to 
other standard models for interacting binaries (Tjemkes et 
al., 1986, Orosz & Bailyn 1997). Briefly, we model a binary 
system in which the primary is a compact invisible object 
and the secondary a ’’normal” star. We assume that the sec¬ 
ondary fills its Roche-lobe, is in synchronous rotation and 
has a circular orbit. A grid consisting of a series of quadrilat¬ 
erals of approximately equal area is set up over the Roche 
surface. The grid is defined using a polar coordinate sys¬ 
tem with Ne rings equally spaced across the star, from the 
Li point to the back of the secondary star (see Orosz & 
Bailyn 1997 for details). For each ring of constant angle 9 
the number of azimuth points is chosen to try to keep the 
area roughly constant. We typically use Ae=40, which cor¬ 
responds to 2048 surface elements. 

The parameters that determine the geometry of the sys¬ 
tem are the binary inclination i, the binary mass ratio q 
(=M 2 /Mi, where Mi and M 2 are the masses of the com¬ 
pact object and secondary star respectively) and the Roche- 
lobe filling factor, /. Since we assume the secondary fills 
its Roche-lobe, a safe assumption in LMXBs, f—1. The ra¬ 
dial velocity semi-amplitude, K 2 , determines the velocity 
scale of the system. The velocity, gravity and temperature 
for each element varies across the star due to the shape of 
the Roche-lobe, the Roche-potential and gravity darkening 
respectively. The temperature and gravity for each element 
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Table 1. SXTS with spectroscopic mass ratio measurements 


Object 

Vrot sin i (km/s) 

q{= M 2 /M 1 ) 

Reference 

GRO J0422-I-32 

90^22 

u.liO_Q 071 

Harlaftis et al., (1999) 

A0620-00 

83±5 

0.067±0.01 

Marsh, Robinson Sz Wood (1994) 

GS2000-I-25 

86±8 

0.042±0.012 

Harlaftis, Horne & Filippenko (1996) 

Nova Mus 1994 

106±13 

n 1 9Q+0.0't4 

Casares et al., (1997) 

Gen X-4 

43±6 

0.17±0.06 

Torres et al., (2002) 

XTE J1550-564 
Nova Sco 1994 

90±10 

82.9-94.9 

0.337-0.436 

Orosz et ah, 2002 

S99 

V4641 Sgr 

123±4 

0.667±0.356 

Orosz et ah, 2001 

V404 Gyg 

39.1±1.2 

u.uuu_o 005 

Casares Si Charles (1994) 


are scaled using the observed effective temperature < Teff > 
and gravity < log g >, which is determined from the spectral 
type and class of the secondary star. 

The local gravity is scaled so that the integrated gravity 
over the star surface is given by the observed gravity. The 
temperature of each element is calculated using the well- 
known von Zeipel relation (von Zeipel 1924), T cc where 
P is the gravity darkening exponent, and is scaled so that the 
integrated luminosity over the stellar surface matches the 
observed bolometric luminosity. The limb angle g = cos^, 
is the angle between the local surface normal and the line of 
sight to the observer; g—1 at the center of the star and g=0 
at the limb. Each element of area on the surface of the star 
is assigned a temperature, gravity, projected velocity and 
limb angle g. For a fixed wavelength, the specific intensity 
for a given element of area depends on Teff, logp and g. 


3.1 The synthetic spectra 

To compute the specific intensity at a given wavelength 
we use the multi-purpose state-of-the-art stellar atmosphere 
code PHOENIX, which produces NextGen model atmo¬ 
spheres, the details of which can be found in Hauschildt, 
Baron & Allard (1997) and Hauschildt et ah, (1999). The 
NextGen models are the most detailed models available 
for cool stars and are computed using spherical geometry, 
rather than the usual plane-parallel approximation. Since 
computing NextGen models are heavily time consuming, 
we cannot compute spectra for all effective temperatures, 
gravities and limb angles. Therefore, we precompute spec¬ 
tra for set Teff, logp and g combinations. 

We compute model spectra with solar abundance in the 
wavelength range 6300 to 6800 A in steps of 0.1 A Tes in the 
range 2400 K to 9800 K in steps of 200 K and loggr in the 
range 1.5 to 5.0 in steps of 0.5. The wavelength range above 
was chosen because this region is commonly used to deter¬ 
mine the secondary star’s rotational broadening since it has 
many strong Fe and Ga absorption line features intrinsic to 
late-type stars. For a given Teff and log g, phoenix computes 
the specific intensities for 64 different angles which are cho¬ 
sen by the code based on the structure of the atmosphere. 
For consistency between the models we interpolate the spe¬ 
cific intensities onto a regular grid with g values ranging 
from 0.05 to 1.0. For accuracy we use a spline interpolation 
method (Press et ah, 1992) in logarithmic units. Also, we 
set the specific intensity to be zero for g < 0.05, since at 
such high limb angles the specific intensity is negligible (see 


Orosz & Hauschildt 2000 for details). Finally, to determine 
the specific intensity for a given wavelength we perform a 3- 
D logarithmic interpolation in Teg, logg and g (see section 
3.2). 

For each element of area visible we compute TeS , log g 
and g. We then determine the specific intensity at each spec¬ 
ified wavelength and integrate the visible specific intensity 
values over the Roche surface. We velocity-correct the model 
spectra by first converting the wavelength scale from vacuum 
to air (Allen 1973) and then bin the spectra onto a uniform 
wavelength scale. Finally, given a colour excess we redden 
the spectrum using the reddening law from of Seaton (1979). 
We thus obtain the exact rotationally broadened spectrum 
of the secondary star at a given orbital phase. The model pa¬ 
rameters that determine the width and shape of the model 
absorption lines are g, i, K 2 and /3; g and K 2 determine 
primarily determine the width of the lines, whereas i and P 
mainly determine the shape of the absorption line (also see 
Shahbaz 1998). 


3.2 The accuracy of the interpolation scheme 

Here we discuss the numerical accuracy of the interpolation 
method we employ. For a fixed wavelength we precompute 
the specific intensity for set of (Teg, loggr, g) grid points. 
Given a model combination of (Teg, logp, g) we perform a 
3-D logarithmic interpolation to determine the specific in¬ 
tensity. Orosz & Hauschildt (2000), showed that the spe¬ 
cific intensity values over a given bandwidth are a relatively 
smooth function of temperature and gravity. We therefore 
tried linear and cubic spline interpolation methods. In con¬ 
trast to Orosz & Hauschildt (200), we found that linear 
or cubic spline interpolation was not accurate and robust 
enough, especially towards the limb of the star. Note that 
we interpolate absorption lines strengths at a given wave¬ 
length, whereas Orosz & Hauschildt (2000) interpolate over 
a large bandwidth. The effects of absorption lines that do not 
vary smoothly would naturally be smoothed out when inte¬ 
grated over a large bandwidth. The interpolation method 
we use is based on a modification of Shepards method for 
a set of 3D scattered data points and is described in Renka 
(1988). Basically the method computes a smooth trivariate 
function using an least squares fit to the grid points near the 
interpolant. This function is then used for the interpolation. 

A simple way to test the interpolation scheme is to re¬ 
move one (Teg, log g and g) grid point from the precomputed 
intensity grid and then calculate its interpolated value. This 
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Figure 1. The effects of using a limb-darkening law to compute 
a Roche-lobe model spectrum. The solid line is the high resolu¬ 
tion model spectrum (top) where the limb darkening is computed 
properly. The model parameters are <7=0.42, it'2=215.5 kms“^, 
i=70°, /3=0.08, < Teff >=6400 K and < logg >=4.0). The 
dashed line is the high resolution model spectrum where a lin¬ 
ear limb-darkening law has been used with a limb-darkening co¬ 
efficient of 0.52. The bottom spectrum is the residual spectrum, 
shifted upwards to clarity. 


value can then be compared to the original grid point value 
and gives an indication of the accuracy of the interpolation 
method. We compute a normal spectrum where we used all 
the grid points in the intensity table and a clipped spec¬ 
trum where we removed the grid point at Tefi^dSOO K, 
log <7=4.5 and /r=0.10 from the table. We find the interpola¬ 
tion method to be accurate to <0.5 percent. To determine to 
effects of using a relatively coarse grid for the Roche-lobe, we 
computed spectra using two different grids, a standard grid 
with Ng—40 and a fine grid with A^e=60, which corresponds 
to 2048 and 4194 grid points respectively. We found that 
the two spectra agreed to within 0.8 percent. Therefore, we 
believe that our integration and interpolation schemes are 
accurate to better than 1 percent. 


3.3 The effects of limb-darkening 

We compare two high resolution Roche-lobe model spectra, 
one where the limb-darkened spectrum is computed prop¬ 
erly and the other where the spectrum is computed asuming 
constant limb darkening (of 0.52) across the line. 

In both cases, the same Roche-lobe system parame¬ 
ters are used (g=0.42, R'2=215.5 kms“^, i=70°, orbital 
phase=0.35, /3=0.08, < Tefj >=6400 K and < logg >=4.0). 
In Fig. one can see that the spectra are clearly differ¬ 
ent. For the constant limb-darkened spectrum, the limb- 
darkening is assumed to be the same for all the lines. Also, 
a linear or two-parameter law does not accurately describe 
the real limb-darkening behavior, since the intensity at the 
limb of the star is much less (Orosz & Hauschildt 2000). 
Therefore it is not surprising that the shape and strength of 
the lines are not the same, especially for the blended lines. 


4 APPLICATION TO THE SXTS: X-RAY 
NOVA SCORPII 1994 

Due to the faintness of the secondary star in most of the 
quiescent SXTs, it is usually only possible to obtain inter¬ 
mediate resolution spectroscopy. However, this implies that 
the shape of the absorption lines cannot be used to deter¬ 
mine system parameters (Shahbaz 1998). The only informa¬ 
tion which we can use is the width of the absorption lines, 
which is related to g. To date there are 9 measurements of 
the binary mass ratio for SXTs (see Table 1), which were 
all obtained by measuring the secondary star’s rotational 
broadening and using the standard limb-darkened rotation 
profile. If we are to determine accurate masses for the bi¬ 
nary components in the SXTs, we have to try and reduce all 
known sources of uncertainties. It is possible to eliminate the 
uncertainties in the binary mass ratio obtained by assuming 
the standard rotation profile, by determining the true limb- 
darkened rotationally broadened spectrum of the secondary 
star. In this section, we outline this method, using the SXT 
Nova Sco 1994 as an example. 


4.1 Previous work 


The SXT Nova Sco 1994 has been extensively studied over 
the past years at both X-rays and optical wavelengths [see 
Shahbaz et ah, 1999 (hereafter S99), Orosz & Bailyn 1997 
and references within]. Optical photometry taken during X- 
ray quiescence (Lx < 10“®Lopt) revealed the orbital period 
Porb=2. 62168 days (Orosz & Bailyn 1997). Through model¬ 
ing the quiescent optical light curves, the binary inclination 
was measured to be 70° (Orosz & Bailyn 1997). Greene, 
Bailyn & Orosz (2001) obtained the same result by model¬ 
ing the optical and infrared quiescent light curves simulta¬ 
neously with an upgraded code. Spectroscopic observations 
allowed the radial velocity semi-amplitude to be measured 
A'2=215.5±2.45 km/s, which implies a binary mass function 
of 2.73TO.O9M0 (S99). The spectral type was determined to 
be F3-F6 (Orosz & Bailyn 1997), and later refined to be F6 
(S99). 

Since the companion star fills its Roche-lobe and is in 
synchronous rotation with the binary motion, Vot sin i com¬ 
bined with K 2 gives a direct measurement of q (Horne, Wade 
& Szkody 1986) 


n sin i = L'2 (1 4- g) 


0.49g' 


,2/3 


0.6g2/3 -kln(l-bg)i/3 


km s"^ (1) 


S99 used the standard technique to estimate the rotational 
broadening of the secondary star. Using optical spectra with 
a velocity resolution of 38.6 kms“^, they compared the av¬ 
erage observed spectrum of Nova Sco 1994 with field tem¬ 
plate stars broadened using the standard rotation profile. 
For limb-darkening coefficients of 0.0 and 0.52, they obtained 
5=0.36 and g=0.40 respectively. However, it should be noted 
that the template star they used has an intrinsic rotational 
broadening of ~30 kms“^(see SIMBAD database), which 
implies that the Vot sin i values were underestimated. In 
order to determine Vot sin i for Nova Sco 1994, we repeated 
the analysis described in S99, but using the F6111 template 
star HR5769 (I4ot sin i<5 kms“^; SIMBAD database) and 
the spectral region 6380A-6500A. For limb-darkening coeffi¬ 
cients of 0.0 and 0.52, we obtain Vot sin i values of 89.6T2.8 
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Wavelength (A) 



Wavelength (A) 


Figure 2. The variance-weighted average spectrum of Nova Sco 1994 (solid line) and the best model spectrum overplotted (dashed line) 
using q=0A2, i=70°, 13=0.08, < Tgff >=6336 K and < logg >=4.0. For the Nova Sco 1994 spectrum we also show the error bars on the 
data. The residual spectrum of Nova Sco 1994 after subtracting the scaled model spectrum is shown underneath. The feature marked 
with an asterisk is residual sky line. The spectra have been normalised. All the spectra are in the rest frame of the model. 


kms“^94.8±3.0 kmrespectively, which correspond to a 
mass ratio of 0.385±0.023 and 0.427±0.024 respectively (90 
percent uncertainty). We find that the spectrum of the sec¬ 
ondary star is not veiled (/„ = !.00±0.05). These results are 
consistent with the Viot sin i and veiling values found by Is- 
raelian et ah, (1999) using high resolution spectra. Note that 
the difference between the Kot sin i values for Nova Sco 1994 
obtained using HR2927 (S99) and HR5769 is ^20 kms“^. 


This difference is the rotational velocity of HR2927 and 
agrees well with the value in the SIMBAD database. 

In order to obtain a more accurate determination of 
q, which does not depend on assumption about the rotation 
profile and limb-darkening coefficients, we use the model de¬ 
scribed in section 3 to compute the exact rotationally broad¬ 
ened model spectrum, which can then be directly compared 
with the average observed spectrum of Nova Sco 1994 pub¬ 
lished in S99. Since the Nova Sco 1994 spectra have a res- 
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olution of 38.6 kms“^, we cannot use the shape of the 
absorption lines to determine system parameters such as i 
or (3 (Shahbaz 1998). The only information which we can 
use is the width of the absorption lines, which is related to 

q- 

4.2 Model parameters 

In our model, the temperature distribution across the sec¬ 
ondary star is described by a gravity darkening law. For 
stars with radiative atmospheres the exponent (/3) for the 
gravity darkening law is 0.25 (von Zeipel 1924), whereas 
if the atmosphere of the star is convective, the exponent 
is 0.08 (Lucy 1967). The spectral type of the secondary 
star in Nova Sco 1994 (F 6 ) is on the boundary between 
hot stars with radiative envelopes and cool stars with con¬ 
vective envelopes. Claret (2000) has shown that the value 
for /3 is a function of the star’s effective temperature. For 
Nova Sco 1994 with < Ten >=6336 K, /3=0.08. However, to 
be conservative in our model computations, we use /3=0.08 
and 0.25. 

We fix i to be 70° (Greene et ah, 2001) and assume 
the secondary to be a F 6 iv (S99) with an < >=6336 K 

and < log ( 7 = 4.0 (Gray 1992). For /3=0.08 the temperature 
and gravity across the secondary star ranges from 5120 K to 
6490 K and log ( 7 = 2.8 to 4.1 respectively. The temperature 
distribution for /3=0.25 is much wider and flatter (3230 K 
to 6780 K) and there is no difference in the range of logg. 

For a compromise between speed and resolution, we 
choose the number of grid points to be such that the maxi¬ 
mum radial velocity difference between any two adjacent ele¬ 
ments is less than the instrumental resolution. We choose the 
total number of surface elements to be 2048 (A^e=40; see sec¬ 
tion 3), which with the binary parameters of Nova Sco 1994, 
gives a maximum radial velocity difference of less than 10 
kms“^. 


4.3 The binary mass ratio 

Our analysis was performed in the same way as by S99, 
the only difference being that we compute broadened spec¬ 
tra using our model which depends on q, instead of us¬ 
ing broadened template stars. For a given value for q we 
first computed the model spectrum in the wavelength range 
6300 A to 6500 A. This region was chosen since it is simi¬ 
lar to the region used in the analysis of S99 and has many 
relatively isolated absorption line features. We then broad¬ 
ened the model spectrum using a Gaussian function with a 
width 38.6 km/s to match the instrumental resolution of the 
Nova Sco 1994 data. We determined the velocity shift of the 
average observed spectrum of Nova Sco 1994 with respect 
to the model spectrum using the method of cross-correlation 
(Tonry & Davis 1979). Prior to the cross-correlation all the 
spectra were rebinned onto a logarithmic wavelength scale 
with a pixel size of 14.5 kms“^. 

In order to compare our model spectrum with the aver¬ 
age observed Nova Sco 1994 spectrum, one should compute 
the model spectra in the same way as the Nova Sco 1994 
spectra were obtained and averaged. The Nova Sco 1994 
average spectrum consisted of 25 velocity-corrected spectra 
taken during 2 nights covering orbital phases near 0.35 and 



Figure 3. The values for the optimal subtraction as a function 
q. The stars and crosses show the values for fits using /3=0.08 
and 0.25 respectively with < Tgfj >=6336 K. The circles show 
the values for fits using < Tgfj >=6600 K and /9=0.08. The 
minimum is 618.6 and the curve has been shifted upwards. 

0.85. Ideally one should compute model spectra at the same 
orbital phases and allow for the orbital smearing due to the 
exposure time. However, the resolution of the data is such 
that any information about the shape of the absorption lines 
is smeared out (Shahbaz 1998). Therefore, we compute a sin¬ 
gle model spectrum at phase 0.35 and assume it is the same 
for all orbital phases. 

To allow for the smearing of the absorption lines due to 
the motion of the secondary star during the length of each 
exposure at each orbital phase, we convolved the model spec¬ 
trum with a rectangular function with the appropriate veloc¬ 
ity width which is given by 2 nK 2 cos{2TT(f))/ kms“^. 
We first normalise using a continuum spline fit and then 
average the 25 model spectra. Finally, to compare the av¬ 
erage model spectrum with the average observed spectrum 
of Nova Sco 1994 we used the method described in Marsh 
et al., (1994). We subtract a constant, /„, multiplied by the 
model spectrum from the average observed Nova Sco 1994 
spectrum. The constant represents the fractional contribu¬ 
tion of the template star spectrum to the total flux density. 
Therefore 1-fv is the veiling factor and for no veiling /„ = 1.0. 
The optimum value for /„ was obtained by minimising the 
statistic. We performed the above analysis for q in the range 
0.32 to 0.4‘52 in steps of 0.01 and in the wavelength range 
6310-6490A. The residual sky feature at 6376A was masked. 
We obtain a minimum of 1729.1 (z.(=530) at (7=0.42 In 
Fig. ^we show the best model spectrum with (7=0.42, the av¬ 
erage observed spectrum of Nova Sco 1994 and the residual 
spectrum. The poor fit is primarily due to the featuress at 
6346A and 6456A contributing significantly to the poor y^. 
Masking out these features we obtain a minimum y^ of 635.3 
((/=463, /3=0.08). Fig. ^ shows the y^ plot for different val¬ 
ues of q. For the models with /3=0.08 (circles) the minimum 
y^ is obtained at q=0.422±0.028 with /„=0.88±0.04 (90 
percent uncertainty). Using a gravity darkening exponent 
of /3=0.25 (crosses), we obtain a significantly higher y^; the 
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Figure 4. The template F6iv star HR5769 (solid line) and the NextGen model spectrum (dashed line) computed using Teff=6400 K 
and log 3=3.7 and degraded to match the resolution of the template star spectrum. The residual spectrum is also shown at the bottom. 
The spectra have been normalised. The wavelength position of the atomic Fe lines and the a-elements are also shown. 


minimum is 646.6 at g=0.439±0.028 with /„=0.89±0.04 
(90 percent uncertainty). 

From optical spectroscopy, S99 and Israelian et al., 
(1999) showed that there is no evidence that the light from 
the secondary star is veiled (see also section 4.1). From our 
fits using < Teft >=6336 K and /1=0.08 we obtain /„=0.88, 
i.e. the model predicts deeper absorption lines compared 
to the data. Since our analysis mainly uses Fe absorption 
lines, which are stronger in cooler stars, we expect by us¬ 
ing a hotter secondary star in the model we will obtain a 


higher value for /„, Using < Tes >=6600 K and /3=0.08 
we obtain a minimum of 618.6 with /„ = 1.00, i.e. the 
secondary star contributes all the flux. The minimum 
is at (j'=0.419±0.028 (90 percent uncertainty). A star with 
< Teff >=6600 K would have an F3 spectral type (Gray 
1992), which is consistent with the observed spectral type 
(F3-F6; Orosz & Bailyn 1997 and S99). 

It should be noted that, although the minimum x^ ob¬ 
tained using < Teft >=6600 K is significantly lower than 
the minimum x^ obtained using < Teft >=6336 K, it should 
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Table 2. Summary of model fits. 90 percent uncertainties are 
given 


Nova Sco 1994 

Standard profile 



u=0.00 

(7=0.385±0.023 


u=0.52 

<?=0.427±0.024 

Nova Sco 1994 

Roche profile 



< Tea >=6336 K, /3=0.08 

(7=0.422±0.028 


< Tea >=6600 K, /3=0.08 

<7=0.419±0.028 

Simulation 

(<7=0.42) 

Standard profile 



u=0.00 

<7=0.392±0.006 


u=0.52 

<7=0.434±0.006 


be noted that the shape and position of the minimum 
does not change; the values are only offset by a constant 
amount. Also, for a fixed < Ten > the fits using /3=0.25 give 
higher values for /„ compared to the fits using /3=0.08 (see 
above). Therefore the determination of q and its uncertainty 
does not depend on the choice of < Tea > • If we use an incli¬ 
nation of 68° or 72°, or if we change K 2 by 4 kms“^(the 90 
percent limits on K 2 ', see S99) the derived value for q does 
not change; the minimum x^ increases by less than 1. Note 
that the uncertainties in the temperature/spectral type and 
in /3 are of the same order as that introduced by using the 
standard rotational profile. 

4-3.1 Overabundant a-elements 

The poor minimum x^ we originally obtained before mask¬ 
ing can be explained either in terms of the NextGen models 
failing to reproduce the strength of some of the absorption 
line features used in the optimal subtraction or that the 
Nova Sco 1994 spectrum contains some overabundant ele¬ 
ments. To see if this mismatch is due to the former, we 
compare a NextGen model spectrum with a F6iv template 
star spectrum. The NextGen model spectrum was com¬ 
puted using refr=6400 K and log ^=4.0, values appropriate 
for the template star HR5769. The model parameters were 
/=0.1 (i.e. for a spherical star), 7^2=0.0 kms“^and /3=0. 
From Fig. ^ it can be seen that we have confidence that our 
code matches the absorption line features: the scatter in the 
residual of the lines is less than 2 percent. 

Israelian et ah, (1999) have measured the metal abun¬ 
dances of the secondary star in Nova Sco 1994 and find 
[Fe/H]=0.1±0.2 which is consistent with being solar. They 
also found an overabundance in the Q-elements (O, Ca, S, 
Mg, Si, Ti and N) compared to solar values. Since the spec¬ 
tral region we use in the optimal subtraction is dominated 
by Fe lines the use of synthetic models with solar abun¬ 
dances seems correct. We would then naturally expect the 
few a-element absorption line features present in the spec¬ 
tral region used for our analysis to be overabundant and 
hence contributes significantly to the poor x^- To deter¬ 
mine the qualitative effects of the a-elements abundances, 
we compare a solar abundance spectrum and a spectrum 
where the a-element abundances have been increased, with 
the observed Nova Sco 1994. (Fig. |^. The abundances of the 
latter spectrum were increased by [O/H] = 1.0, [S/H]=0.75, 


[Mg/H]=0.9, [Si/H]=0.9, [Ti/H]=0.9 and [N/H]=0.45 (Is¬ 
raelian et ah, 1999). As Israelian et ah, (1999) did not 
determine the Ga abundance, we also increased the Ga 
abundance by [Ca/H]=1.0 to see if the mismatched fea¬ 
ture at 6456A is due to Ga. Both spectra were computed 
using reff=6400 K and log 5=4.0 and were broadened by 
Kot sin i=95 kms“^using the standard rotation profile 
(u=0.52). Subtracting solar- and enhanced-abundance spec¬ 
tra from the observed Nova Sco 1994 spectrum, we ob¬ 
tain Xv values of 2.95 and 4.15 respectively. The residual 
solar- and enhanced-abundance spectra are shown in Fig. ^ 
The mismatched feature at 6436A contain the a-elements 
Mgii 6346.742 A and Si 6415.48 A, whereas the feature at 
6456A contains Gal 6455.598 A and Oi 6455.977 A but note 
that Fell 6456.38 A is also present. 

While our spectra seem to show that some of the mis¬ 
matched features are due to the a-elements being being over¬ 
abundant, there are also a few features that the non-solar 
spectrum predicts to be overenhanced. To some extent this 
explains why the x2 for the enhanced spectrum is higher 
than the xS for the solar spectrum. Given the poor resolu¬ 
tion of the spectra, it is difficult to determine with certainty 
if the mismatched features are entirely due to the a-elements 
being overabundant. Only a more detailed abundance analy¬ 
sis using high resolution spectra can prove this convincingly. 

5 DISCUSSION 

5.1 The mass ratio and black hole mass 

Beer & Podsiadlowski (2002) determined q by fitting the 
optical lights curves of Nova Sco 1994 in a self consistent 
manner, without any assumptions about the distance or in¬ 
terstellar reddening. They obtained g=0.256±0.039 which 
is significantly lower than the values obtained by S99 (see 
also section 4.1) and Greene et ah, (2001). They argue that 
the high mass ratio obtained by Greene et ah, (2001) is 
a consequence of neglecting the accretion disc when mod¬ 
eling the optical and infrared light curve. Furthermore, in 
order to explain the disagreement in q, they suggested that 
the secondary star’s rotational broadening measured by S99 
was overestimated. [Note that the revised determination of 
q (see section 4.1) does not affect the argument presented 
by Beer & Podsiadlowski (2002)]. They comment that since 
the spectra in S99 were taken at quadratures near orbital 
phases 0.25 and 0.75 (orbital phase 0.0 is defined as supe¬ 
rior conjunction of the secondary star) the effects of orbital 
smearing would be the greatest. The orbital smearing would 
add a linear amount of 3.5 kms“^to the rotational broad¬ 
ening. As S99 did not take into account the effects of orbital 
smearing, the rotational broadening value obtained by S99 
would be overestimated. 

Firstly, we would like to point out that the effects of 
orbital smearing would be the least because the velocity of 
the secondary star is least at the quadratures. Secondly, be¬ 
cause the velocity width of the rectangular smearing profile 
is much smaller than the rotational broadening, the effects 
of orbital smearing is very small. To examine the effects of 
orbital smearing on the rotational broadening we performed 
a simple simulation where we computed a Gaussian func¬ 
tion with FWHM=100 kms“^, appropriate for the observed 


© 0000 RAS, MNRAS 000, 000-000 


Determining the spectroscopic mass ratio in interacting binaries 9 



6320 6340 6360 6380 6400 

Wavelength (A) 



6400 6420 6440 6460 6480 

Wavelength (A) 


Figure 5. A comparison between a solar and non-solar LTE spectrum. The LTE synthetic spectrum were computed with < >=6400 K 

and log 3=4.0. The solid line shows the LTE spectrum with solar abundance. The dashed lines shows the LTE spectrum where the ex- 
elements have abundances of [O/H]=1.0, [S/H]=0.75, [Mg/H]=0.9, [Si/H]=0.9, [Ti/H]=0.9, [N/H]=0.45 (Israelian et ah, 1999) and 
[Ca/H]=1.0. The LTE model spectra have been broadened by 95 kms“^. The residual spectrum after subtracting the solar (middle) 
and enhanced-abundance (bottom) model spectrum from the observed Nova Sco 1994 spectrum are also shown. 

Kot sin i (see Royer et al., 2002 for a rule of thumb calibra- tematic uncertainties in modeling the light curves (e.g. the 

tion) and then smeared it with a rectangular profile with binary inclination is correlated with Teff) it is difficult to ex- 

a width of 7 kms“^. We find that the width of the Cans- plain the difference in q obtained by Beer & Podsiadlowski 

sian function before and after the smearing is the same to (2002). 

within 0.2 percent, therefore, the effects of orbital smearing 
are negligible. Also, it should be noted that Beer & Podsi¬ 
adlowski (2002) obtained Tefr^OlSO K which is significantly Using our value for the binary mass ratio, we can de- 

different to what we obtain (see section 4.3). Given the sys- termine the mass of the binary components from the mass 

function equation 
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f{M) = 


PKi _ Ml sin® i 
2^ ^ (l + g)2 




( 2 ) 


With f(M)=2.73±0.15 Mg (S99), (j=0.419±0.028 (see sec¬ 
tion 4.3) and i=70.2±1.6° (Greene et al., 2001), we obtain 
Mi= 6.59±0.45 Mg and M2=2.76±0.33 Mg. The 90 percent 
nncertainites were obtained using a Monte Carlo procedure. 
Although the the absolute value for the binary mass ratio 
obtained by Greene et ah, (2001) and in this paper are dif¬ 
ferent, the uncertainties in the mass ratio are such that the 
black hole mass we obtain is consistent with that obtained 
by Greene et ah, (2001). 


5.2 Roche-lobe versus standard rotation profile 

In S99 we determined q for Nova Sco 1994 by comparing 
the average observed spectrum of Nova Sco 1994 with the 
spectrum of a template star which had been convolved with 
the standard rotation profile. The standard rotation profile 
requires a limb-darkening coefficient which we assume to be 
the value appropriate for the continuum i.e. we assume that 
the radiation in the lines is the same as for the continuum. 
Since the absorption lines in early type stars will have core 
limb-darkening coefficients much less than the values appro¬ 
priate for the continuum (Collins & Truax 1995), in order 
to be conservative we also used a limb-darkening coefficient 
of zero. Using a non-rotating template star with line limb- 
darkening coefficients of zero and 0.52 (note that 0.52 is 
actually the continuum value), we obtained mass ratio’s of 
g=0.385 and 0.427 respectively (see section 4.1). Using the 
Roche-lobe model we obtained q=0.419 (see section 4.3). As 
one can see, the Roche-model gives a mass ratio that lies 
in between the values obtained using limb-darkening coeffi¬ 
cients of zero and 0.52 (continuum) and the standard rota¬ 
tion profile. 

To prove this, we simulate a spectrum and then try 
to determine Kot sin i and hence q in the same way as 
we did in S99 (see section 4.3) and using the X-ray binary 
model (section 4.1). We hrst perform the same analysis as in 
S99 to determine Mot sin i, i.e. we compare broadened ver¬ 
sions of the non-rotating template star HR5769 with the 
simulated average spectrum, using the standard rotation 
profile with a continuum limb-darkening coefficients of 0.0 
and 0.52. The high quality (S/N=500) model spectrum was 
computed using <7=0.42, 7^2=215.5 kms“^, i=70°, /1=0.08, 
< Teff >=6400 K and < log g >=4.0. For limb-darkening co¬ 
efficients of 0.0 and 0.52 we obtain mass ratio’s of q=0.392 
and 0.434 respectively (90 percent uncertainty of 0.006). Us¬ 
ing the binary model, we hnd that we recover the exact 
model q value. These simulations conhrm that, for the case 
of an F-type secondary star, the standard rotation profile 
with zero and continuum value for the line limb-darkening 
coefficient gives a value for q that brackets the value found 
using the X-ray binary model. 

The limb-darkening coefficient varies across the absorp¬ 
tion having the continuum value in the wings to a reduced 
value in the line core, where the exact value for the line 
limb-darkening coefficient depends on the atomic species of 
the line. Since most of the absorption in the rotationally 
broadened line comes from the line core, the value for the 
limb-darkening in the line core will dominate the shape of 
the line (Collins & Truax 1995). The fact that we obtain 


a better estimate for q using a limb-darkening coefficient 
slightly less than the continuum value suggests that the av¬ 
erage line limb-darkening coefficient for the strongest ab¬ 
sorption lines used in the analysis (i.e. Fe) is close to the 
continuum value. It is clear that only detailed line profile 
calculations on individual lines can determine the true ab¬ 
sorption line limb-darkening coefficient. 


6 CONCLUSIONS 

We present a model that predicts the spectrum of the sec¬ 
ondary in an interacting binary. The model uses synthetic 
NextGen model spectra which are incorporated into the 
secondary star’s Roche geometry. As a result, we determine 
the exact rotationally broadened spectrum of the secondary 
which does not depend on assumptions about the rotation 
profile and limb-darkening coefficients. 

As an example, we determine the mass ratio for the 
SXT Nova Sco 1994. We use our model to compute the ro¬ 
tationally broadened model spectrum, which we compare 
directly with the observed intermediate resolution spectrum 
of Nova Sco 1994. In order to match the width of the ab¬ 
sorption lines and the lack of veiling in the quiescent optical 
spectrum, our model requires q=0.419±0.028 (90 percent 
uncertainty), and a secondary star with < Teff >=6600 K 
(an F3 star) which is consistent with the observed spec¬ 
tral type. For the case of an F-type secondary star, we Hnd 
that the standard rotation profile with zero and continuum 
value for the line limb-darkening coefficient gives a value for 
q that brackets the value found using the full geometrical 
treatment. 
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